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Abstract 

We use molecular dynamics simulations to probe the rotational dynamics of the SPC/E model of 
water for a range of temperatures down to 200 K, 13 K above to the mode coupling temperature. 
We find that rotational dynamics is spatially heterogeneous, i.e., there are clusters of molecules 
that rotate significantly more than the average for a given time interval, and we study the size and 
the temporal behavior of these clusters. We find that the position of a rotational heterogeneity 
is strongly correlated with the position of a translational heterogeneity, and that the fraction of 
molecules belonging to both kinds of heterogeneities increases with decreasing temperature. We 
further find that although the two types of heterogeneities are not identical, they are related to 
the same physical picture. 

PACS numbers: 61.20.Ja, 61.20.Gy 
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Experiments and computer simulations have shown that dynamics in supercooled liq- 
uids is spatially heterogeneous, i.e^one can identify transient regions with relaxation times 
different by orders of magnitude 111. Simulations have shown that the most mobile parti- 
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cles tend to form clusters M, . Different theoretical approaches have been developed to 
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understand spatially heterogeneous dynamics j^|5|. In particular Adam and Gibbs (AG) 
0, 0] postulate the existence of cooperatively rearranging regions (CRR) whose molecules 
change configuration independently of the rest of the system. Molecular dynamics (MD) 
simulations p, have verified many of the predictions of the AG theory, which has also 
been tested in other systems j^. 

While there have been numerous studies of the heterogeneous nature of the translational 
degrees of freedom (TDOF) in water, there are few studies of the heterogeneous nature of the 
rotational degrees of freedom (RDOF). Here we systematically study the rotational dynamics 
of water and compare with the translational dynamics. Previous work for other systems 
sug gests that translationally mobile molecules may have enhanced rotational mobilitv 0, 
Ul, and the characteristic times for the RDOF are smaller than for the TDOF To 
this end we perform MD simulations of a system of = 1728 water molecules interacting 
with the extended simple point charge potential (SPC/E) for a range of temperatures 
from 350 K down to 200 at the fixed density of 1 g/cm^; for each temperature we run two 
independent trajectories to improve the statistics. Simulation details are given in jl^ . 

We quantify the rotation of a molecule using the normalized polarization vector Pi{t). In a 
time interval [t,t + 6t], pi{t) spans an angle 69 = cos^^ {Pi{t) " Pi{t + ^t))- We define a vector 
60i{t) such that = 66 and its direction is given by Pi{t) x pi[t + 6t) . Thus the vector 

ifii^t) = /q 6(pi{t')dt' allows us to define a trajectory in a (y9-space representing the rotational 
motion of molecule i. One can then associate a rotational mean square displacement (RMSD) 
given by 

i 

and a rotational diffusion coefficient 

1 ^ 

= }^ E( I^^W - ^^(O)l') • (2) 

i=l 

The vector 0i{t) is not bounded to the unit sphere, since otherwise Eq. Q would give 
Dr = 0. Equations ((T)) and ^ were applied in [l5| to study the rotational motion of a linear 
molecular system. 



Figures Ufa) andmb) show the RMSD and the temperature dependence of Dr, respec- 
tively. Similar to Ref. we observe: (i) the RMSD shows three different regimes: a 
ballistic regime, where < ip"^ >oc t"^, a plateau or cage regime, where molecules find them- 
selves trapped in the cage formed by their neighboring molecules, and finally a diffusive 
regime where < ip"^ >oc t, these three regimes are analogous to those observed in studies of 
translational dynamics of supercooled liquids, (ii) Dr increases with T with a non-Arrhenius 
behavior. Note the oscillations in Fig. ^a) present for times of the order of 10~^ ps. These 
oscillations, not present for linear molecules correspond to the lihration (hindered ro- 
tation) regime and occur at the same time as the oscillations observed in the rotational 
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correlation function of water in jl6l ]. 

We introduce the rotational counterpart of the self part of the time-dependent van Hove 
distribution function, Gs((y9,t) [l^, 

1 ^ 

t) ^ -Y.^5{\m - - ^)) , (3) 

i=l 

where (■ ■ ■) represents average over configurations. With this formalism we recover the usual 
interpretation for 47r(y9^Gs((y9, t) as the probability of having a molecule at time t with angular 
displacement (p. In other words, that in the abstract <y9-space, a molecule has moved to a 
distance (p> from its position at t = 0. For long times the diffusion equation for (pi{t) holds, 
and Gs{ip,t) is a Gaussian distribution 
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exp [-3^y2{v\t))] . (4) 



The deviations of Gs{'p,t) from Go{ip,t) can be quantified by the non-Gaussian parameter 

Q 

FigureEl^a) shows 02 (^) for different temperatures. a2{t) shows a clear maximum at t = t*(T) 
which corresponds to the beginning of the diffusive regime for the RDOF [Fig. ^a)]. We 
note that there is a small maximum at t ~ 10"^ ps: this is a consequence of the librational 
motion as shown in Fig. d^a) Figure |2l^c) shows Gsifjt) and Go{ip,t) for T = 200 K 
and t = t*{200K) ^ 1.05 ns. As in [sj we find that G^i^^t*) and Go{<p,t*) intersect for 
large ip at ip*, and Gs^ipyt*) shows a large tail where the fitted Gaussian underestimate the 
angular motion of the molecules. Molecules with ip > ip* can be considered with an angular 
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displacement higher than expected; this fraction / = Anip'^Gs^ip, t) is found to be ~ 13%, 
showing a weak T- dependence. 



In analogy to 



2l|, we define the rotational mobility of a molecule at a given time 



to as the maximum angular displacement in the interval [to,to + of th^ oxygen atom 

^iit, At) = max{|^i(t + to) - 0^{to)\ , to < t < to + At} . (6) 

We focus our attention on the most rotationally mobile molecules and explore the possibility 
that there exist clusters also among this category of molecules. To facilitate comparison with 



the study 



of translational heterogeneities (TH), we select the 7% of the most rotationally 



mobile molecules j2^ and define a cluster at time to over an observation time At as those 
molecules whose nearest neighbor oxygen-oxygen (0-0) distance at time to is less than 
0.315 nm (first minimum of 0-0 radial distribution function). We find that the rotational 
dynamics is spatially heterogeneous, since these molecules form clusters, which we will call 
rotational heterogeneities (RH). 

Let us now address the question of how these RH depend on the observation time. Let 
{n{At)) be the average number of molecules in a cluster a time At. We define the weight 
average cluster size as 

To eliminate the contribution of random clusters, we normalize (n(At))^„ by {nr)w, i-e., the 
weight average cluster size obtained by chosing randomly 7% of the molecules. Figure El 
shows {n{At))w/{nr)w as a function of At for different T. In the same manner as the 
translational case, the maximum in {n^At))^/ {nr)w occurs at the end of the cage regime of 
the RMSD, indicating that the cage breaking of the RDOF is highly correlated with the cage 
breaking of TDOF. The RH become larger as T decreases. FigurelH^a) shows {n{At))w/ {nr)w 
for RH and TH at the corresponding tmax, the time at which the corresponding weight 
average cluster size is largest. We find that on average TH are larger than RH, and that RH 
reach their maximum size before the TH, Fig. Efb), i.e. tmax for RH is smaller than tmax for 
TH. 

It is natural to ask if RH and TH are formed by the same molecules. We address this 
question by simultaneously analyzing the properties of RH and TH. Figure Efa) is a typical 
snapshot of the system, showing both TH and RH. The clusters together form a larger 
entity characterizing the dynamical heterogeneities; the molecules belonging to both kinds 



of clusters act as the "backbone" of such an entity. We find that the fraction of molecules 
simultaneously belonging to both clusters increases with decreasing temperature. 

To compare the structure of the TH and RH, we evaluate the radial distribution function 
(RDF) of oxygen atoms within each cluster, and between the two clusters. Figure El^b) shows 
the RDF for TH, for RH, and for molecules which belong to both TH and RH. We see that 
there is a strong tendency for mobile molecules (of either type) to be neighbors. The RDF's 
are qualitatively similar to the bulk RDF, gbuik{i^), with maxima at r ^ 0.28 nm and r ^ 
0.45 nm (i.e., molecules are nearest or next-nearest neighbors); however, the amplitudes of 
the first peak are strongly enhanced compared to bulk water. In order to examine deviations 
from the bulk we normalize the RDF's by gbuik{f) [inset of Fig.j^b)]. All of the RDF's display 
maxima at 0.32 nm, indicating that oxygens in TH and RH have an enhanced tendency (with 
respect to the bulk) to be in the first interstitial shells of each other and, therefore, have more 
than four nearest neighbors. Molecules with five or more neighbors have bifurcated bonds 
and represent "defects" in the tetrahedral network characterizing water [22, 24|. Therefore, 
our results suggest that the combined TH and RH in water (i) are a consequence of the 
defects of the HB network, and (ii) are primarily composed by molecules located at the 
defects of the HB network: evidence of the relation between the structure and dynamics 
in water. Thus the physical picture needed to describe rotational heterogeneities resembles 
that needed for describing translational heterogeneities. 
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FIG. 1: (a) RMSD for a range of temperatures from 200 to 350 K. (b) Rotational diffusivity as a 
function of T^^. 
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FIG. 2: (a) Non-Gaussian parameter 02 (i) for the range of temperatures indicated, (b) T- 
dependence of t*{T), the time at which the maximum of 02 occurs. For comparison, we also 
show the the value of t* obtained from the study of the TDOF in P]. (c) Self part of the van 
Hove distribution function, Gs{f,t), for T = 200 K and t*(200 K) ~ 1.05 ps, compared with the 
Gaussian approximation, Go{(p,t), obtained using {ip'^{t*)) also at T = 200 K. 
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FIG. 3: Weight average cluster size for T = 200 K — 350 K for molecules belonging to the 7% most 
rotationally mobile molecules. The values are normalized to the random cluster contribution. 
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FIG. 4: (a) Weight average cluster size as a function of temperature for TH and RH. Both quantities 
are calculated at the corresponding tmax- Note that the TH are larger on average than the RH. 
(b) The time t^ax at which the maximum of the weight average cluster size occurs. 
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FIG. 5: (a) A snapshot of one cluster at T = 200 K. The red molecules belong to TH, blue to RH 
and the green ones belong to both clusters, (b) RDF between the sets of translationally mobile 
and rotationally mobile molecules for T = 200 K. The inset shows the corresponding ratios of these 
functions to g{r)buik- For clarity qt-t is shifted by 5 and qr-r by 10 on the vertical axis; g^_rp 
by 3 and g%^_R by 6 in the inset. 
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